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A first search is reported for a standard model Higgs boson (H) that is produced 
through vector boson fusion and decays to a bottom-quark pair. Two data samples, 
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for a Higgs boson mass of 125 GeV. 
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1 Introduction 

In the standard model (SM) ID-EI^ the electroweak symmetry breaking is achieved by a mech¬ 
anism 12-01 that provides mass to the electroweak gauge bosons, while leaving the photon 
massless. The mechanism predicts the existence of a scalar Higgs boson (H), and its observa¬ 
tion was one of the main goals of the CERN LHC program. A boson with mass near 125 GeV 
was recently discovered by both the ATLAS |Zl and CMS 13 H] collaborations, with properties 
that are compatible with those of a SM Higgs boson llOllHl . 

At the LHC, a SM Higgs boson can be produced through a variety of mechanisms. The expected 
production cross sections 1121 as a function of the Higgs boson mass are such that, in the mass 
range considered in this study, the vector boson fusion (VBL) process pp —> qqH has the second 
largest production cross section following gluon fusion (GL). Lurthermore, for a SM Higgs 
boson with a mass nin < 135 GeV, the expected dominant decay mode is to a b-quark pair (bb). 

Thus far, the search for H —> bb has been carried out in the associated production process 
involving a W or a Z boson (VH production) at the Tevatron IT^ and at the LHC IIT2 [151 , as 
well as in association with a top quark pair at the LHC IIT61 - IT81 , without reaching the necessary 
sensitivity to observe the Higgs boson in this decay channel. It is therefore important to exploit 
other production modes, such as VBL, to provide in the bb decay channel further information 
on the nature and properties of the Higgs boson. 

The prominent feature of the VBL process qqH ^ qqbb is the presence of four energetic jets 
in the final state. Two jets are expected to originate from a light-quark pair (u or d), which are 
typically two valence quarks from each of the colliding protons scattered away from the beam 
line in the VBL process. These "VBL-tagging" jets are expected to be roughly in the forward 
and backward directions relative to the beam direction. Two additional jets are expected from 
the Higgs boson decay to a bb pair in more central regions of the detector. Another impor¬ 
tant property of the signal events is that, being produced through an electroweak process, no 
quantum chromod}mamics (QCD) color is exchanged at leading order in the production. As 
a result, in the most probable color evolution of these events, the VBL-tagging jets connect to 
the proton remnants in the forward and backward beam line directions, while the two b-quark 
jets cormect to each other as decay products of the color neutral Higgs boson. Consequently 
very little additional QCD radiation and hadronic activity is expected in the space outside the 
color-cormected regions, in particular in the whole rapidity interval (rapidity gap) between the 
two VBL-tagging jets, with the exception of the Higgs boson decay products. 

The dominant background to this search is from QCD production of multijet events. Other 
backgrounds arise from: (i) hadronic decays of Z or W bosons produced in association with 
additional jets, (ii) hadronic decays of top quark pairs, and (iii) hadronic decays of singly pro¬ 
duced top quarks. The contribution of the Higgs boson in GL processes with two or more 
associated jets is included in the expected signal yield. 

The search is performed on selected four-jet events that are characterized by the response of a 
multivariate discriminant trained to separate signal events from background without making 
use of kinematic information on the two b-jet candidates. Subsequently, the invariant mass 
distribution of two b jets is analyzed in each category in the search for a signal "bump" above 
the smooth contribution from the SM background. This is the first search of this kind, and the 
only search for the SM Higgs boson in all-jet final states at the LHC. A search for a SM Higgs 
boson in the all-hadronic final state has been previously reported by the CDL experiment Jill. 

This paper is organized as follows: Section|^highlights the features of the CMS detector needed 
to perform this analysis. Section details the production of simulated samples used to study 
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3 Simulated samples 


the signal and main backgrounds, and Section [^presents the employed triggers. Event recon¬ 
struction and selection are described in Sections and respectively The unique features of 
the analysis are discussed in Section which includes the improvement of the resolution in 
jet transverse momentum (pj) by regression techniques, discrimination between quark- and 
gluon-originated jets, and soft QCD activity An important validation of the analysis strategy 
is the observation of the Z —> bb decay, which is presented in Section The search for a SM 
Tiiggs boson is discussed in Section|^and the associated systematic uncertainties are presented 
in Section]^ The final results are discussed in Sectionj^and combined with previous searches 
in the same charmel in Section We summarize in Section 


2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. A silicon pixel and strip tracker, a lead tungstate crystal 
electromagnetic calorimeter, and a brass and scintillator hadron calorimeter are located within 
the solenoidal field. Muons are measured in gas-ionization detectors embedded in the steel 
flux-return yoke of the solenoid. Forward calorimetry (pseudorapidity 3 < \r]\ < 5) comple¬ 
ments the coverage provided by the barrel (|7/| < 1.3) and end cap (1.3 < \ri\ < 3) detectors. 
The first level (LI) of the CMS trigger system, composed of specialized processors, uses infor¬ 
mation from the calorimeters and muon detectors to select the most interesting events in a time 
interval of less than 4 ps. The high-level trigger (HLT) processor farm decreases the event rate 
from about 100 kHz to less than 1 kHz, before data storage. A more detailed description of the 
CMS apparatus and the main kinematic variables used in the analysis can be found in Ref. Il20ll . 


3 Simulated samples 

Samples of simulated Monte Carlo (MC) signal and background events are used to guide the 
analysis optimization and to estimate signal yields. Several event generators are used to pro¬ 
duce the MC events. The samples of VBF and GF signal processes are generated using the next- 
to-leading order perturbative QCD program POWHEG 1.0 HZTl . interfaced to TAUOLA 2.7 1221 
and PYTHIA 6.4.26 ||23l for the hadronization process and modeling of the underlying event 
(UE). The most recent PYTHIA 6 Z2* tune is derived from the Z1 tune If24l , which uses the 
CTEQ5L parton distribution functions (PDF), whereas Z2* adopts CTEQ6L ESil . The signal 
samples are generated using only H —i bb decays, for five mass hypotheses: nin = 115, 120, 
125,130, and 135 GeV. 

Background samples of QGD multijet, Z +jets, W+jets, and tf events are simulated using leading- 
order (LO) MadGraph 5.1.3.2 Il26l interfaced with PYTHIA. The single top quark background 
samples are produced using POWHEG, interfaced with TAUOLA and PYTHIA. The default set 
of PDF used with POWHEG samples is GTIO Il27l . while the LO GTEQ6L1 set Il25l is used for 
other samples. The production cross sections for W+jets and Z +jets are rescaled to next-to- 
next-to-leading-order (NNLO) cross sections calculated using the FEWZ 3.1 program Il28l - l30l . 
The tt and single top quark samples are also rescaled to their cross sections based on NNLO 
calculations 113111521 . 

To accurately simulate the LHG luminosity conditions during data taking, additional simulated 
pp interactions overlapping in the same or neighboring bunch crossings of the main interac¬ 
tion, denoted as pileup, are added to the simulated events with a multiplicity distribution that 
matches the one in the data. 
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4 Triggers 

The data used for this analysis were collected using two different trigger strategies that result in 
two different data samples for analysis. First, a set of dedicated trigger event selection (paths) 
was specifically designed and deployed for the VBF qqH —> qqbb signal search, both for the 
LI trigger and the HLT, and operated during the full 2012 data taking. Then, a more general 
trigger was employed for the larger part of the 2012 data taking that targeted VBF signatures 
in general. The first (nominal) set of triggers collected the larger fraction of the signal events, 
while the second trigger supplemented the search with events that failed the stringent nominal- 
trigger requirements. The integrated luminosity collected with the first set of triggers was 
19.8 fb^^, while for the second trigger it was 18.3 fb^^. 

While the first dedicated trigger paths collected data within the standard CMS streams, the sec¬ 
ond general-purpose VBF trigger path ran in parallel with data streams that were reconstructed 
later, in 2013, during the LHC upgrade. 

4.1 Dedicated signal trigger 

The LI paths require the presence of at least three jets with pj above decreasing thresholds 
Pj \ pj \ pj'^ that were adjusted according to the instantaneous luminosity [pj^ = 64-68 GeV, 
p^'^ = 44-48GeV, pj'^ = 24-32GeV ]. Among the three jets, one and only one of the two 

leading jets [with px > Pj \ Px^^ ] can be in the forward region with pseudorapidity 2.6 < IpI < 
5.2, while the other two jets are required to be central (|p | < 2.6). 

The HLT paths are seeded by the LI paths described above, and require the presence of four jets 
with Px above thresholds that were again adjusted according to the instantaneous luminosity, 
Px > 75-82, 55-65, 35M8, and 20-35 GeV, respectively. Two complementary HLT paths have 
been employed that make use, respectively, of (i) only calorimeter-based jets (Galojets) and (ii) 
particle-flow jets (PFJets, see Section]^. At least one of the four selected jets must further fulfill 
minimum b-tagging requirements, evaluated using HLT regional tracking around the jets, and 
using the "track counting high-efficiency" (TGHE) or the "combined secondary vertex" (GSV) 
algorithms |j33l, alternatively for the first and second paths. Events are accepted if they satisfy 
either of the two paths. 

Among the four leading jets, the light-quark (qq) VBF-taggmg jet pair is assigned in one of 
two ways: (i) the pair with the smallest HLT b-tagging output values (b-tag-sorted qq) or (ii) 
the pair with the maximum pseudorapidity difference (p-sorted qq). Both pairs are required 
to exceed variable minimum thresholds on |Apqq| of 2.2-2.5, and of 200-240 GeV on the dijet 
invariant mass niqq, depending on the instantaneous luminosity. 

To evaluate trigger efficiencies, a prescaled control path is used, requiring one PFJet with px > 
80 GeV. To match the efficiency in data, the simulated trigger efficiency must be corrected with 
a scale factor of order 0.75 that is parametrized a function of the highest jet b-tag output in 
the event and the invariant mass of the quark-jet candidates. With this procedure the weak 
dependence of the trigger efficiencies on the invariant mass of the two b jets is also taken into 
account. 

4.2 General-purpose VBF trigger 

The LI paths for the general-purpose VBF trigger require that the scalar px sum of the hadronic 
activity in the event exceeds 175 or 200 GeV, depending on the instantaneous luminosity. 

The HLT path is seeded by the LI path described above, and requires the presence of at least 
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5 Event reconstruction 


two Calojets with pj > 35 GeV. Out of all the possible jet pairs in the event, with one jet lying 
at positive and the other at negative p, the pair with the highest invariant mass is selected as 
the most probable VBF-tagging jet pair. The corresponding invariant mass and absolute 
pseudorapidity difference are required to be larger than 700 GeV and 3.5, respectively. 

The efficiency of the general-purpose VBF trigger is measured in a similar way as for the ded¬ 
icated triggers, using a prescaled path (requiring two PFJets with average pj > 80 GeV). To 
match the efficiency in data, the simulated trigger efficiency must be corrected with a scale fac¬ 
tor of order 0.8 that is expressed as a function of the invariant mass and the pseudorapidity 
difference of the two offline quark-jet candidates. 


5 Event reconstruction 

The offline analysis uses reconstructed charged-particle tracks and candidates from the particle- 
flow (PF) algorithm Il34l - t36ll . In the PF event reconstruction all stable particles in the event, i.e. 
electrons, muons, photons, and charged and neutral hadrons, are reconstructed as PF candi¬ 
dates using information from all GMS subdetectors to obtain an optimal determination of their 
direction, energy, and type. The PF candidates are then used to reconstruct the jets and missing 
transverse energy. 

Jets are reconstructed by clustering PF candidates with the anti-kj algorithm Il37ll38ll with a dis¬ 
tance parameter of 0.5. Reconstructed jets require a small additional energy correction, mostly 
due to thresholds on reconstructed tracks and clusters in the PF algorithm and various recon¬ 
struction inefficiencies (391. Jet identification criteria are also applied to reject misreconstructed 
jets resulting from detector noise, as well as jets heavily contaminated with pileup energy (clus¬ 
tering of energy deposits not associated with a parton from the primary pp interaction) (401. 
The efficiency of the jet identification criteria is greater than 99%, with the rejection of 90% of 
background pileup jets with pj ~ 50 GeV. 

The identification of jets that originate from the hadronization of b quarks is done with the GSV 
b tagger |33l , also implemented for the HLT paths, as described in SectionThe GSV algo¬ 
rithm combines the information from track impact parameters and secondary vertices identi¬ 
fied within a given jet, and provides a continuous discriminator output. 

Events are required to have at least four reconstructed jets. All the jets found in an event are 
ordered according to their pj, and the four leading ones are considered as the most probable 
b jet and VBF-tagging jet candidates. The distinction between the two jet types is done by the 
means of a multivariate discriminant that, in addition to the b-tag values and the b-tag order¬ 
ing, takes into account the rj values and the p ordering. In the VBF H —)■ bb signal simulation it 
is found that the b jets have higher b-tag values and are more central in p than the VBF-tagging 
jets. A boosted decision tree (BDT), implemented with the TMVA package (4TI . is trained on 
simulated signal events using the discriminating variables previously described and its output 
is used as a b-jet likelihood score; out of the four leading jets the two with the highest score are 
identified as the b jets, while the other two are identified as the VBF-tagging jets. With the use 
of the multivariate b-jet assignment the signal efficiency is increased by ^10% compared to the 
interpretation based on GSV output only. 
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6 Event selection 

The offline event selection is based upon the b-jet and VBF-tagging jet assignment described in 
Section]^ and is adjusted to the two different trigger sets presented in Section]^ In what follows 
the selected events are divided into two sets referred to as set A and set B. These selections are 
summarized in Table [H 

Events selected in set A are required to have been selected by the dedicated VBF qqH —> qqbb 
trigger and to have at least four PF jets with > 80,70,50,40 GeV and \ t]\ <4.5. Moreover, 

at least two of these jets must satisfy the loose CSV working point requirement (CSVL) Il33l . 
The VBF topology is ensured by requiring m^q > 250 GeV and |A?^qq| > 2.5, where qq denotes 
the pair of the most probable VBF-tagging jets. Finally, in order to suppress further the back¬ 
ground, the azimuthal angle difference A^b between the two b-jet candidates must be less than 
2.0 radians. Figurel^shows the normalized distributions of |A;yqq| (left) and A<^b (right) for the 
sum of all simulated backgrounds, and the VBF and GF Higgs boson production. 

Events in set B are first required to not belong to set A (to avoid double counting). Then, 
they must have passed the generic VBF topological trigger and have at least four PF jets with 
Pt > 30 GeV and | p | < 4.5. In addition, the scalar pj sum of the two leading jets must be greater 
than 160 GeV. In order to enrich the sample in b jets, there must be at least one jet satisfying 
the medium CSV working point requirement (CSVM) 1331 and one jet satisfying the CSVL. The 
VBF topology is ensured by requiring rriqq, > 700GeV, and \Arjqq\, > 3.5, where 

qq denotes the pair of the most probable VBF jets and jj denotes the jet pair with the highest 
invariant mass (as in the trigger logic described in Section]^. Finally the azimuthal angle A(^b 
between the two b-jet candidates must be less than 2.0 radians. 

After all the selection requirements, 2.3% of the simulated VBF signal events (for = 125 GeV) 
end up in set A, and 0.8% end up in set B. The fraction of events in set A that also satisfy the 
requirements of set B (except for the set A veto) amounts to 39%. The set B selection recovers 
signal events presenting pronounced VBF jets, with lower pj but larger pseudorapidity open¬ 
ing and invariant mass. 


Table 1: Summary of selection requirements for the two analyses. 



Set A 

SetB 

Trigger 

Dedicated VBF qqH ^ qqbb 

General-purpose VBF trigger 



p™ > 30 GeV 

Jets Pt 

> 80,70,50,40 GeV 

Px -h Px > 160 GeV 

Jets IpI 

<4.5 

<4.5 

btag 

At least 2 CSVL jets 

At least 1 CSVM and 1 CSVL jets 

A^b 

<2.0 radians 

<2.0 radians 


mqq> 250 GeV 

mqq, > 700 GeV 

VBF topology 

Apqq > 2.5 

lApqql, lApf^l > 3.5 

Veto 

None 

Events that belong to set A 
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7 Signal properties 



• Data (set A) 
-VBF H(125} 
..GFH(125) 

[] Background 
[^MC stat. unc. 


(a) 


CMS 


19.8fb ' (8TeV) 



Figure 1: (a) Normalized distribution in absolute pseudorapidity difference between the two 
VBF-jet candidates (|A)^qq|). (b) Normalized distribution of the azimuthal difference between 
the two b-jet candidates (A(/)bb). The selection corresponds to set A, data are shown by the 
points, and the sum of all simulated backgrounds is by the filled histograms. The VBF Higgs 
boson signal is displayed by a solid line, and the GF Higgs boson signal is shown by a dashed 
line. The panels at the bottom show the fractional difference between data and background 
simulation, with the shaded band representing the statistical uncertainties in the MC samples. 


7 Signal properties 

The analysis described in this paper relies on certain characteristic properties of the studied 
final state, which provide a significant improvement of the overall sensitivity. First, the reso¬ 
lution of the invariant mass of the two b jets is improved by applying multivariate regression 
techniques. Then, the jet composition properties are used to form a discriminant that separates 
jets originating from light quarks or gluons. Third, soft QCD activity outside the jets is quan¬ 
tified and used as a discriminant between QCD processes with strong color flow and the VBF 
signal without color flow. 


7.1 Jet transverse-momentum regression 

The bb mass resolution is improved by using a regression technique similar to those used in the 
search for a Higgs boson produced in association with a weak vector boson and decaying to 
bb mu. A refined calibration is carried out for individual b jets, beyond the default jet energy 
corrections, that takes into account the jet composition properties and targets semileptonic b 
decays that lead to a substantial mismeasurement of the jet pj due to the presence of an es¬ 
caping neutrino. For this purpose a regression BDT is trained on simulated signal events with 
inputs including information about the jet properties and structure. The target of the regres¬ 
sion is the Pj of the associated particle-level jet, clustered from all stable particles (with lifetime 
CT > 1 cm). The inputs include: (i) the jet pj, rj, and mass; (ii) the jet energy fractions carried 
by neutral hadrons and photons IIMM36ll : (iii) the mass and the uncertainty in the decay length 
associated with the secondary vertex, when present; (iv) the event missing transverse energy 
and its azimuthal direction relative to the jet; (v) the total number of jet constituents; (vi) the pj 
of the soft-lepton candidate inside the jet, when present, and its pj component perpendicular 
to the jet axis; (vii) the pj of the leading track in the jet; and (viii) the average pj density of the 
event in y-cp space (FastJet p method HU). 

The additional energy correction of b jets leads to an improvement of the jet pj resolution, 
which in turn improves the dijet invariant mass resolution by approximately 17% in the phase 



















7.2 Discrimination between quark- and giuon-originated jets 


7 


space of the offline event selections. Figure shows the reconstructed dijet invariant mass of 
the b-jet candidates (wtbb) before and after the regression for simulated events passing set A 
selections. The measured distribution of the regressed mbb in set A is shown in Fig. 

The validation of the regression technique in data is done with samples of Z —events 
with one or two b-tagged jets. When the jets are corrected by the regression procedure, the px 
balance distribution, between the Z boson, reconstructed from the leptons, and the b-tagged jet 
or dijet system is improved to be better centered at zero and narrower than when the regression 
correction is not applied. In both cases the distributions for data and the simulated samples are 
in good agreement after the regression correction is applied. 


CMS Simulation 



Figure 2: Simulated invariant mass distribution of the two b-jet candidates before and after the 
jet Px regression, for VBF signal events. The generated Higgs boson signal mass is 125 GeV and 
the event selection corresponds to set A. By FWHM we denote the width of the distribution at 
the middle of its maximum height. 

7.2 Discrimination between quark- and giuon-originated jets 

To further identify whether the jet pair with the smallest b-tag values among the four leading 
jets is likely to originate from hadronization of a light (u,d,s-type) quark, as expected for signal 
VBF jets, or from gluons, as is more probable for jets produced in QCD processes, a quark-gluon 
discriminant Il43ti45l is applied to the VBF candidate jets. 

The discriminant exploits differences in the showering and fragmentation of gluons and quarks, 
making use of the following internal jet composition observables based on the PF jet con¬ 
stituents: (i) the major root-mean square (RMS) of the distribution of jet constituents in the 
rj-cp plane Il45l . (ii) the minor RMS of the distribution of jet constituents in the rj-cp plane Il45l . 
(iii) the jet asymmetry pull |46l, (iv) the jet particle multiplicity, and (v) the maximum energy 
fraction carried by a jet constituent. The pull and RMS variables are calculated by weighting 
each jet constituent by its px squared 1451 . 

The five variables above are used as inputs to a likelihood estimated with gluon and quark jets 
from simulated QCD dijet events using the TMVA package. To improve the separation power, 
all variables are corrected for pileup effects as a function of the FastJet p density. Figure 
shows the normalized distribution of the quark-gluon likelihood (QGL) I45l for the first VBF- 











8 
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Regressed m (GeV) 

bb 


• Data (set A) 

□ QCD (x 1.68} 

[]Z+jets 

□tt 

[]Single t 
[]W+jets 
-lOxVBF H(125) 
-- lOxGF H(125) 
|]MC stat. unc. 


Figure 3: Distribution in invariant mass of the two b-jet candidates, after the jet px regression, 
for the events of set A. Data are shown by the points, while the simulated backgrounds are 
stacked. The LO QCD cross section is multiplied by a factor 1.68 so that the total number of 
background events equals the number of events in the data, while the VBF and GF Higgs boson 
signal cross sections are multiplied by a factor 10 for better visibility The last bin is the sum 
of all the events beyond the range of the x axis (overflow). The panel at the bottom shows the 
fractional difference between the data and the background simulation, with the shaded band 
representing the statistical uncertainties in the MC samples. 


jet candidate (the jet that is ranked third in the b-tag score; see Section]^, for background and 
signal events. As expected, VBF signal events, dominated by quark jets, have a pronounced 
peak at likelihood ~0, while the background and GF events are enriched in gluon jets, and 
have a very different QGL distribution. The QGL distribution of all four jets is used as input to 
the signal vs background discriminant (Section |9.1[ ). 

7.3 Soft QCD activity 

To measure the additional hadronic activity between the VBF-tagging jets, excluding the more 
centrally produced Higgs boson decay products, only reconstructed charged tracks are used. 
This is done to measure the hadronic activity associated with the primary vertex (PV), defined 
as the reconstructed vertex with the largest sum of squared transverse momenta of tracks used 
to reconstruct it. 


A collection of "additional tracks" is assembled using reconstructed tracks that (i) satisfy the 
high purity quality requirements defined in Ref. Il47l and pj > 300 MeV; (ii) are not associated 
with any of the four leading PF jets in the event; (iii) have a minimum longitudinal impact pa¬ 
rameter, |dz(PV) |, with respect to the main PV, rather than to other pileup interaction vertices; 
(iv) satisfy |d 2 (PV)| < 2mm and |(iz(PV)| < 3(7'z(PV) with respect to the PV, where UzjPV) 
is the uncertainty in dz{PY); and (v) are not in the region between the two best b-tagged jets. 
This is defined as an ellipse in the f]-(p plane, centered on the midpoint between the two jets, 
with major axis of length AR(bb) -|- 1, where AR(bb) = (A//bb)^ + (A^b)^/ oriented along 

the direction cormecting the two b jets, and with minor axis of length 1. 

The additional tracks are then clustered into "soft Trackjets" using the anti-kj clustering algo¬ 
rithm with a distance parameter of 0.5. The use of Trackjets represents a clean and validated 
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19.8fb ' (8TeV) 



• Data (set A) 
-VBF H(125} 
--GFH(125} 

[] Background 
|]mc stat. unc. 


Figure 4: Normalized distribution in quark-gluon likelihood discriminant of the first light-jet 
candidate. Quark jets are expected to have low likelihood values (closer to 0), while gluon 
jets are expected to have higher ones (closer to 1). The selection corresponds to set A, data 
are shown by the points, and the sum of all simulated backgrounds is shown by the filled his¬ 
togram. The VBF Higgs boson signal is displayed by a solid line, and the GF Higgs boson signal 
is shown by a dashed line. The panel at the bottom shows the fractional difference between the 
data and the background simulation, with the shaded band representing the statistical uncer¬ 
tainties in the MC samples. 

method 1481 to reconstruct the hadronization of partons with very low energies down to a few 
GeV [Ml; an extensive study of the soft Trackjet activity can be found in Refs. l4^l44ll . 

The discriminating variable, that encapsulates the differences between the signal and the 

QGD background, is the scalar pj sum of the soft Trackjets with pj > 1 GeV, and is shown in 

Kg-@ 


8 Extraction of the Z boson signal 

The Z +jets background process, with the Z boson decaying to a b-quark pair, provides a vali¬ 
dation of the analysis strategy used for the Higgs boson search. The extraction of the Z boson 
signal demonstrates the ability to observe a relatively wide hadronic resonance on top of a 
smooth QGD background. Also, if such a signal can be seen, it can serve for in situ confirma¬ 
tion of the scale and resolution of the invariant mass of the two b jets. Recently, the observation 
of a Z —> bb signal was reported by the ATLAS Gollaboration 1501 in the Z +1 jet final state, 
and similar techniques are applied here. The overall strategy has two parts. First, events are 
selected from set A, with the additional requirement to have at least one GSVM jet. It should 
be noted that it is important to extract the Z boson signal in the same four-jet phase space in 
which the Higgs boson search is performed. Then, a multivariate discriminant is trained to 
separate the Z +jets process from the QGD multijet production, using variables that are only 
weakly correlated to rntb- According to the value of the discriminant the events are divided 
into three categories, ranging from a signal-depleted control category, to a signal-enriched one. 
Finally, a simultaneous fit of the signal and the QGD background f«bb shape is performed in all 
three categories. The subsequent sections give details of the outlined procedure. 
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8 Extraction of the Z boson signal 



• Data (set A) 
-VBF H(125} 
--GFH(125} 

[] Background 
|]mc stat. unc. 


Figure 5: Normalized distribution of the scalar px sum of Trackjets that are associated with the 
soft QCD activity The selection corresponds to set A, data are shown by the points, and 

the sum of all simulated backgrounds is shown by the filled histogram. The VBF Higgs boson 
signal is displayed by a solid line, and the GF Higgs boson signal is shown by a dashed line. 
The panel at the bottom shows the fractional difference between the data and the background 
simulation, with the shaded band representing the statistical uncertainties in the MC samples. 


8.1 Z boson signal vs. background discrimination 

As discussed above, the selection of events is based on set A, with the additional requirement 
of having at least one CSVM jet; the tightening of the b-tagging condition was found to im¬ 
prove the expected sensitivity. A Fisher discriminant (FD) I14TI is implemented with the TMVA 
package and trained to discriminate between the Z +jets signal and the background. For this 
purpose, seven variables are used: (i) the absolute rj difference | Apqq | of the VBF jets; (ii) the ab¬ 
solute Tj of the b-jet system |pbb |i (hi) the CSV value of the jet with highest CSV value (with best 
b tag); and (iv)-(vii) the QGL values of the four leading jets. Due to the very small correlations 
between the variables, the FD performs almost as well as more advanced, nonlinear discrim¬ 
inators. Figure shows the normalized distribution of the discriminant, where the output of 
the Z +jets signal is compared to the background. 


8.2 Fit of the dijet invariant mass spectrum 

The selected events are divided into three categories, based on the FD output. Table [^summa¬ 
rizes the event categories and corresponding yields. 

Table 2: Definition of the event categories for the Z boson signal extraction and corresponding 
yields in the ?«bb interval [60,170] GeV. 



Category 1 

Category 2 

Category 3 


FD < -0.02 

-0.02 < FD < 0.02 

FD > 0.02 

Data 

659873 

374797 

342931 

Z +jets 

1374 

1467 

2783 

tt 

2124 

1821 

2327 

Single t 

657 

569 

812 
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CMS 


19.8fb ' (8TeV) 



• Data (set A) 
— Z+jets 
[] Background 
|]MC stat. unc. 


Figure 6: Normalized distribution in Z boson Fisher discriminant. Data are shown by the 
points, and the sum of all simulated backgrounds is shown by the filled histogram. The Z 
+jets signal is displayed with solid line. The panel at the bottom shows the fractional differ¬ 
ence between the data and the background simulation, with the shaded band representing the 
statistical uncertainties in the MC samples. 

Besides its discrimination power, the FD has minimal correlation with the invariant mass of 
the two b jets. This means that the ?«bb spectrum from QCD processes is independent of the 
category In practice, however, there is a small residual dependence, of up to 3%, which is 
corrected with a linear transfer function of wjbb that is taken from data. The extraction of the 
Z boson signal is done with a simultaneous fit in all three categories. Eq. Q describes the fit 
model: 

fii^hh) = FzM',zZj(?«bb;kjES/kjER) + N;,tT;(?«bb) + Ni^QCD B8(?«bb;F)/ (1) 

where the subscript i denotes the category, N, z is the expected yield for the Z boson signal; and 
Fz/ M,qcd are free parameters for the signal strength and the QCD event yield. The shape of the 
top quark background T; (?«bb) is taken from the simulation (sum of the tt and single top quark 
contributions), and the expected yield N;q is allowed to vary in the fit by 20%. The Z +jets 
signal shape Z;(wJbbi^jES/^jER) is taken from the simulation and is parametrized as a Crystal 
Ball function ll^ (Gaussian core with power-law tail) on top of a combinatorial background 
modeled by a polynomial. The position and the width of the Gaussian core are allowed to vary 
by the factors kjEs and /cjer, respectively, which quantify any mismatch of the jet energy scale 
and resolution between data and simulation and are constrained by the dedicated validation 
measurements of the regressed jet energy scale and resolution. 

Finally, the QGD background shape in each category is described by a common, eighth-order 
polynomial B8(mbb;p)/ whose parameters p are determined by the fit, and a multiplicative 
linear transfer function ]<C; (mbb) that accounts for the small background shape differences be¬ 
tween the categories. The choice of the polynomial is based on an extensive bias study de¬ 
scribed in Section Allowing for 20% uncertainty on the Z boson signal efficiency, the si¬ 
multaneous binned maximum-likelihood fit yields a signal strength of pz = which 

corresponds to an observed (expected) significance of 3.6 a (3.3 cr). The fitted values of /cjes and 
kjER are 1.01 ± 0.02 and 1.02 ± 0.10, respectively Figure [^shows the fitted distributions and the 
background-subtracted ones. 
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9 Search for a Higgs boson 


The extraction of the Z boson signal in this way validates the Higgs boson search method used 
in this paper by finding a known dijet resonance in a similar mass range. In addition, the 
simulated ?nbb scale and resolution are consistent with the data, based on the best-fit values of 
the kj^s arid /cjer nuisance parameters, which serve to constrain the corresponding uncertainties 
in the Higgs boson signal extraction. 


X10^ 19.8fb '(8TeV) x10^ 19.8 fb'(8 TeV) 



(GeV) 


Figure 7: Invariant mass distribution of the two b-jet candidates for the Z boson signal in 
the three event categories that are based on the Z boson Fisher discriminant output, start¬ 
ing from the most backgroundlike (upper left) and ending at the most signal-like (bottom). 
Data are shown by the points. The solid line is the sum of the postfit background and signal 
shapes, while the dashed line is the background component alone. The bottom panel shows 
the background-subtracted distribution, overlaid with the fitted signal, and with the Icr and 2cr 
background uncertainty bands. The measured (simulated) parameters of the Gaussian core of 
the signal shape in category 3 are 97.7 (96.6) GeV and 9.3 (9.1) GeV for the mean and the sigma, 
respectively. 


9 Search for a Higgs boson 

The search for a Higgs boson follows closely the methodology applied for the extraction of 
the Z boson signal (Section |8.2| . Namely, a multivariate discriminant is employed (Section |9.1| 
to divide the events into seven categories that are subsequently fit simultaneously with 
templates (Section [9^. 
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9.1 Higgs boson signal vs. background discrimination 

In order to separate the overwhelmingly large QCD background from the Higgs boson signal, 
all discriminating features have to be used in an optimal way. This is best achieved by using a 
multivariate discriminant, which in this case is a BDT implemented with the TMVA package. 
The variables used as an input to the BDT are chosen such that they are very weakly correlated 
with the d}mamics of the bb system, in particular with m^b/ and are grouped into five distinct 
groups: (i) the dynamics of the VBF-jet system, expressed by Ar^qq, Acpqq, and niqq, (ii) the b-jet 
content of the event, expressed by the CSV output for the two best b-tagged jets; (iii) the jet 
flavor of the event QGL for all four jets; (iv) the soft activity, quantified by the scalar pj sum 
of the additional soft Trackjets with px > 1 GeV, and the number of soft Trackjets 
with px > 2 GeV; and (v) the angular dynamics of the production mechanism, expressed by 
the cosine of the angle between the qq and bb planes in the center-of-mass frame of the four 
leading jets cos0*q^,^. In practice, two BDTs are trained with the same input variables using 
the selections corresponding to the two sets of events. This distinction is necessary because the 
properties of the selected events are significantly different between the two selections (set A 
and set B). Figure|^shows the output of the BDT for the two sets of events. 


CMS 19.8fb-'(8TeV) 



BDT output 


• Data (set A) 
□QCD (xl.65) 
□Z+jets 
□ti 

□Singlet 
□w+jets 
-lOxVBF H(125) 
-- lOxGF H{125) 
□ mC stat. unc. 


(a) 


CMS l8.3fb'(8TeV) 



BDT output 


• Data (set B) 
□QCD (X 1.80) 
□Z+jets 
□tf 

□ Singlet 
□w+jets 
-10xVBFH(125) 
-■ 10xGFH(125) 

□ mC stat. unc. 

(b) 


Figure 8: Distribution of the BDT output for the events of set A (a) and set B (b). Data are shown 
by the points, while the simulated backgrounds are stacked. The LO QGD cross sections are 
scaled such that the total number of background events equals the number of events in data, 
while the VBF and GF Higgs boson signal yields are multiplied by a factor of 10 for better 
visibility. The panels at the bottom show the fractional difference between the data and the 
background simulation, with the shaded band representing the statistical uncertainties of the 
MG samples. 


9.2 Fit of the dijet invariant mass spectrum 

Taking into account the expected sensitivity of the analysis and the available number of MG 
events (necessary to build the various templates), seven categories are defined, according 
to the BDT output: four for set A and three for set B. The boundaries of the categories and 
the respective event yields are summarized in Table In an mbb interval of twice the width 
of the Gaussian core of the signal distribution (wjh = 125 GeV), the signal-over-background 
ratio reaches 1.7% in the most sensitive category (category 4). It should be noted that both the 
VBF and GF contributions are added to the Higgs boson signal, with the fraction of the latter 
ranging from ~50% in category 1 to ~7% in category 4. 

The analysis relies on the assumption that the QGD mbb spectrum shape is the same in all 
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9 Search for a Higgs boson 


Table 3: Definition of the event categories and corresponding yields in the m\,\, interval 
[80,200] GeV, for the data and the MC expectation. The BDT output boundary values refer 
to the distributions shown in Fig. 


BDT boundary values 

set A 

set B 

Cat. 1 

Cat. 2 

Cat. 3 

Cat. 4 

Cat. 5 

Cat. 6 

Cat. 7 

O 

O 

1 

kO 

O 

1 

0.0-0.7 

0.7-0.84 

0.84-1.0 

-0.1-0.4 

0.4-0.8 

0.8-1.0 

Data 

546121 

321039 

32740 

10874 

203865 

108279 

15151 

Z +jets 

2038 

1584 

198 

71 

435 

280 

45 

W+jets 

282 

135 

4 

<1 

225 

92 

17 

tf 

2818 

839 

45 

14 

342 

169 

21 

Single t 

960 

633 

64 

25 

194 

159 

30 

VBF mH(125) 

53 

140 

58 

57 

33 

57 

31 

GF mH(125) 

53 

51 

8 

5 

9 

10 

2 


BDT categories of the same set of events. In reality, a small correction is needed to account for 
residual differences between the ?«bb spectrum in category 1 vs. categories 2,3 and 4, and in 
category 5 vs. categories 6 and 7. The correction factor (transfer function) is a linear function 
of wJbb hr set A and a quadratic one in set B (because a stronger dependence is observed in 
set B between wjbb arid the multivariate discriminant). With the introduction of the transfer 
functions, the fit model for the Higgs boson signal is given by Eq. (|^: 

fi{^hh) = Eh H,(mbb;^JES/%R) + Z;(/«bb;^JES/%R) 

+M',t T';(mbb;^JES/%ER) + Ni^QCD E(?nbb; pset), 

where the subscript i denotes the category and }iUf M,qcd are free parameters for the signal 
strength and the QCD event yield. N, Z/ and are the expected yields for the Higgs 

boson signal, the Z +jets, and the top quark background respectively. The shape of the top 
quark background T, (mbb;kjES/^jER) is taken from the simulation (sum of the tt and single 
top quark contributions) and is described by a broad Gaussian. The Z/W+jets background 
Z;(mbb; kjES/ ^jer) and the Higgs boson signal H;(mbb; kjES/ ^jer) shapes are taken from the sim¬ 
ulation and are parametrized as a Grystal Ball function (Gaussian core with power-law tail) on 
top of a polynomial. The position and the width of the Gaussian core of the MG templates (sig¬ 
nal and background) are allowed to vary by the free factors /cjes and Ejer, respectively, which 
quantify any mismatch of the jet energy scale and resolution between data and simulation. 
Finally, the QGD shape is described by a polynomial B(mbb;pset), common within the cate¬ 
gories of each set, and a multiplicative transfer function fC, (?7rbb) per category, accounting for 
the shape differences between the categories. The parameters of the polynomial, p^eu and those 
of the transfer functions, are determined by the fit, which is performed simultaneously in all 
categories in each set. For set A, the polynomial is of fifth order, while for set B it is of fourth 
order. 

The choice of the QGD background shapes and event category transfer function parametriza- 
tions are fully based on data, and have been performed among classes of functions, e.g. poly¬ 
nomials, exponential, power laws and their combinations, with a minimum number of degrees 
of freedom suited to fit the data in all categories. Each function considered is used to generate 
different MG pseudo-data sets, and each data set is fitted using the different functional models. 
A potential bias on the signal estimation is computed for each pair of possible functions used 
to generate and fit to the pseudo-data sets. The background model chosen yields a maximum 
potential bias on the fitted signal strength of less than six times the statistical uncertainty on 
the background. Hence the systematic uncertainty associated with the background shape can 
be neglected. 
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10 Systematic uncertainties 

Table [^summarizes the sources of uncertainty related to both the background and to the signal 
processes. The leading uncertainty comes from the QCD background description: both the pa¬ 
rameters of its shape and the overall normalization in each category are allowed to float freely, 
being determined by the simultaneous fit to the data. The resulting covariance matrix is used 
to compute the uncertainty. For the smaller background contributions from the Z/W+jets and 
top quark production, the mbb shapes are taken from the simulation, while their correspond¬ 
ing yields are allowed to float in the fit with a 30% log-normal constraint centered on the SM 
expectations. 

The experimental uncertainties on the jet energy scale (JES) and jet energy resolution (JER) 
affect the signal acceptance and the shape of the multivariate discriminant output, and are in¬ 
cluded as nuisance parameters. The effect of the JES and JER uncertainties on the m\,\, shape 
is taken into account in the fit function. By varying the JES and JER by their measured uncer¬ 
tainties |[39l| , the impact of the signal yield per analysis category is estimated. These variations 
affect the acceptance by up to 10%, while the peak position of the mbb shape is shifted by 2%, 
and the width by 10%. 

Additional uncertainties are assigned to the flavor tagging of the jets. The CSV and QGL dis¬ 
criminant outputs are shifted according to the observed agreement between data and simu¬ 
lation and the effect on the signal acceptance is estimated to range from 3% to 10% for the 
former, and from 1% to 3% for the latter. The impact of the CSV shift is more significant, both 
because it is used for the event selection, and because the multivariate discriminant depends 
more strongly on the b tagging of the jets. The shift of QGL only affects the shape of the dis¬ 
criminant, and thus the distribution of signal events in the analysis categories. 

The trigger uncertainty is estimated by propagating the uncertainty in the data vs. MG simu¬ 
lation scale factor for the efficiency. This is achieved by convolving the two-dimensional effi¬ 
ciency scale factor with the signal distribution. As a result, the uncertainty in the signal yield 
ranges from 1% to 6% for the VBF process, and from 5% to 20% for the GF. 

Theoretical uncertainties affect the signal simulation. First, the uncertainty due to PDFs and 
strong coupling constant as variation is computed to be 2.8% (VBF) and 7.5% (GF) l52l . A resid¬ 
ual uncertainty from these sources is estimated for the particular kinematical phase space of the 
search: following the PDF4LHG prescription 11531154| the PDF and as uncertainty ranges from 
2% to 5%, while the renormalization and factorization scale variations in the signal simulation 
induce an uncertainty of 1% to 5% in the analysis categories, on top of a global cross section un¬ 
certainty of 0.2% (VBF) and 7.7-8.1% (GF). Finally, the variation of the UE and parton-shower 
(PS) model (using PYTHIA 8.1 |[55l instead of the default pythia 6) affects the signal acceptance 
by 2% to 7% (VBF) and by 10% to 45% (GF). 

Lastly, an uncertainty of 2.6% is assigned to the total integrated luminosity measurement Il5^ . 


11 Results 

The J^bb distributions in data, for all categories, are fitted simultaneously with the paramet¬ 
ric functions described in Section 9.2 under two different hypotheses: background only and 
background plus a Higgs boson signal. The fit is a binned likelihood fit incorporating the sys¬ 
tematic uncertainties discussed in Section [Tol as nuisance parameters. Due to the smallness 


of the GF contribution in the most signal-sensitive categories we do not attempt to fit inde¬ 
pendently the VBF and the GF signal strengths. The fits in sets A and B are shown in Figs. [^ 
and respectively. The limits on the signal strength are computed with the asymptotic GLg 
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12 Combination with other CMS Higgs boson to b-quarks searches 


Table 4: Sources of systematic uncertainty and their impact on the shape and normalization of 
the background and signal processes. 


Background uncertainties 

QGD shape parameters 

determined by the fit 

QGD bkg. normalization 

determined by the fit 

Top quark bkg. normalization 

30% 


Z/W+jets bkg. normalization 

30% 


Uncertainties affecting the signal 

VBF signal GF signal 

JES (signal shape) 

2% 


JER (signal shape) 

10% 


Integrated luminosity 

2.6% 


Branching fraction (H —)■ bb) 

2.4%M.3% 

JES (acceptance) 

6%-10% 

4%-12% 

JER (acceptance) 

1%M% 

l%-9% 

b-jet tagging 

3%-9% 

3%-10% 

Quark/gluon-jet tagging 

l%-3% 

l%-3% 

Trigger 

l%-6% 

5%-20% 

Theory uncertainties 

VBF signal GF signal 

UE&PS 

2%-7% 

10%M5% 

Scale variation (global) 

0.2% 7.7%-8.1% 

Scale variation (categories) 

l%-5% 

l%-5% 

PDF (global) 

2.8% 

7.5% 

PDF (categories) 

1.5%-3% 

3.5%-5% 


method IISTh^ . Figure 11 shows the observed (expected) 95% confidence level (CL) limit on 
the total VBF plus GF signal strength, as a function of the Higgs boson mass, which ranges from 
5.0 (2.2) at wjh = 115 GeV to 5.8 (3.7) at ?«h = 135 GeV, together with the expected limits in the 
presence of a SM Higgs boson with a mass of 125 GeV. For the 125 GeV Higgs boson signal the 
observed (expected) significance is 2.2 (0.8) standard deviations, and the fitted signal strength 
is = a /ugM = 2.8^1'^. The measured signal strength is compatible with the SM Higgs boson 
prediction ^ = 1 at the 8% level. 


12 Combination with other CMS Higgs boson to b-quarks searches 


The GMS experiment has also performed searches for the Higgs boson decaying to bottom 
quarks, where the Higgs boson is produced in association with a vector boson ^T41 (VH), or 
with a top quark pair lH^llTl (tfH). The VH results have been recently updated and combined 
with tfH llTT]| . Here we combine those results with the ones from the VBF production search 
described in this paper. Event selection overlaps between different analyses have been checked 
and are either empty by construction or have negligible effects on the combination. The com¬ 
bination methodology is based on the likelihood ratio test statistics employed in Section 
and takes into account correlations among sources of systematic uncertainty. Gare is taken to 
understand the behavior of the parameters that are correlated between analyses, in terms of 
the fitted parameter values and uncertainties. 
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Table lists the 95% GL expected and observed upper limits and the best-fit signal strength 
values from the individual charmels and from the combined fit. For nin = 125 GeV the com¬ 
bination yields an H —)■ bb signal strength = l.Od^g^^ with a significance of 2.6 standard 
deviations. 
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Figure 9: Fit of the invariant mass of the two b-jet candidates for the Higgs boson signal (mu = 
125 GeV) in the four event categories of set A. Data are shown by the points. The solid line 
is the sum of the postfit background and signal shapes, the dashed line is the background 
component, and the dashed-dotted line is the QCD component alone. The bottom panel shows 
the background-subtracted distribution, overlaid with the fitted signal, and with the Ic and 2cr 
background uncertainty bands. 
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13 Summary 
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Figure 10: Fit of the invariant mass of the two b-jet candidates for the Higgs boson signal 
{mil = 125 GeV) in the three event categories of set B. Data are shown with markers. The solid 
line is the sum of the postfit background and signal shapes, the dashed line is the background 
component, and the dashed-dotted line is the QCD component alone. The bottom panel shows 
the background-subtracted distribution, overlaid with the fitted signal, and with the \a and 2a 
background uncertainty bands. 


13 Summary 

A search has been carried out for the SM Higgs boson produced in vector boson fusion and 
decaying to bb with two data samples of pp collisions at \/s = 8 TeV collected with the CMS 
detector at the LHC corresponding to integrated luminosities of 19.8 fb^^ and 18.3 fb~^. Upper 
limits, at the 95% confidence level, on the production cross section times the H —?> bb branching 
fraction, relative to expectations for a SM Higgs boson, are extracted for a Higgs boson in the 
mass range 115-135 GeV. In this range, the expected upper limits in the absence of a signal 
vary between a factor of 2.2 to 3.7 of the SM prediction, while the observed upper limits vary 
from 5.0 to 5.8. For a Higgs boson mass of 125 GeV, the observed and expected significance 
is, respectively, 2.2 and 0.8 standard deviations, and the fitted signal strength is }i = a /ctsm = 
2.8^^'^. This is the first search of this kind, and the only search for the SM Higgs boson in all-jet 
final states, at the LHC. 

The combination of the results obtained in this search with other CMS H —?► bb searches in the 
VH and ttH production modes, yields a H —> bb signal strength pi = l.OS+g^^ with a signal 
significance of 2.6 standard deviations for wth = 125 GeV that is consistent with the SM. 
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Figure 11: Expected and observed 95% confidence level limits on the signal cross section in 
units of the SM expected cross section, as a function of the Higgs boson mass, including all 
event categories. The limits expected in the presence of a SM Higgs boson with a mass of 
125 GeV are indicated by the dotted curve. 


Table 5: Observed and expected 95% CL limits, best fit values on the signal strength parameter 
ji = a /ctsm and signal significances for = 125 GeV, for each H —> bb channel and their 
combination. 


H^bb 

Channel 

Best fit (68% CL) 
Observed 

Upper limits (95% CL) 
Observed Expected 

Signal significance 
Observed Expected 

VH 

ttH 

VBF 

0.89 ± 0.43 

0.7 ±1.8 

1.68 

4.1 

5.5 

0.85 

3.5 

2.5 

2.08 

0.37 

2.20 

2.52 

0.58 

0.83 

Combined 

i.uo o.42 

1.77 

0.78 

2.56 

2.70 
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